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Abstract — The recent Auger results suggest that although 
coincidences of arrival directions with 'nearby' AGN and HiRes 
discovery of the GZK cut-off indicate protons, the measured 
longitudinal propagation characteristics indicate heavy nuclei, 
if the conventional interaction model is correct. Something has 
to change! Our own view is that it is possible that the AGN 
-implied proton identification is not correct and that the extra 
galactic particles are, in fact, mainly 'heavies', in which case the 
interaction problem goes away. However, here we assume that the 
particles ARE protons and examine the possible consequences. 
Parameters discussed include the interaction mean-free-path, 
inelasticities and 'exotic' possibilities. 

I. Introduction 

Even a very brief review of connections between the Cosmic 
Ray (CR) physics and physics of high energy interactions has 
to start from the origin of the High Energy Physics (HEP) 
which has taken place high in the atmosphere where the 
multiparticle production processes involving CR particles have 
been seen in nuclear emulsions. 

The second important conjunction point is associated with 
the first and simplest statistical model describing statistically 
these phenomena proposed by Fermi which was falsified also 
in the emulsion balloon experiments showing the asymmetry 
known today as the jet structure. 

Then, for some time, it was believed that the pre-ISR 
(Intersecting Storage Rings) HEP has reached an asymptotia 
and nothing unexpected appear when the energy rise. The 
cross sections were looking stable, according to the theory, 
as well as average transverse momenta, which were close to 
the Hagedorn limit of the highest ever possible temperature 
[1]. However, even in this pacific time of the end of '60 
there were suppositions, experimental indications (see, e.g., 
[2]) from cosmic ray physics where the energies were orders 
of magnitude higher than available for accelerator HEP that 
the Extensive Air Showers (EAS) do not look like they should, 
if nothing unexpected appears there. This suspicions were 
undoubtedly confirmed by ISR. The cross sections started 
to increase, the transverse momenta distributions rise tails 
showing - just unexpected. The ISR revolution initiated, to 
some extent, by the CR discoveries, drove to a fall of the 
thermodynamic, Hagedorn interaction picture (not entirely of 
course) and an ascent of the parton/quark hadronization related 
to the name of Feynman [3]. 

The Feynman scaling worked well and the quark hadroniza- 
tion picture, beautiful and simple satisfying everyone - almost 
everyone. The debate concerning the Feynman scaling (in the 
forward region) initiated again in the cosmic ray community 



was resolved (with the certainty possible in colliding exper- 
iment when the very forward region is unseen) by the next 
machine, the Super Proton Synchrotron (SPS). This happened 
about the a quarter of century ago, just around the time when 
there was the European Cosmic Ray Symposium in Kosice for 
the first time. 

These four facts are of course a subjective choice, but, I 
think, that they ought to find considerable places in any list 
of the interaction points between the CR and HE Physics. 

The first remark I would like to underline in this introduc- 
tion is: the long time has been passed since last event on the 
list. 

The second point is related to attempt to answer the ques- 
tion, why this is so. Is it true that CR physics really has 
no discovery potential for contemporary HEP? The wishful 
answer is, of course, "not!" 

In this paper we would like not only to diagnose the 
situation and propose the treatment, but also to show an 
example of the effect of such therapy. We will show the 
prediction to be tested in the future (not so very far future, 
indeed) related to discoveries recently announced by the Pierre 
Auger Observatory [4] and HiRes [5]. 

II. CORSIKA 

The last decade of the previous century, was very successful 
in the particle physics, mainly in the theory. In the field of 
cosmic rays the attempt was made to answer one of the most 
important questions: about the nature of 'the knee' in the CR 
energy spectrum. 

The question was stated precisely in Kernforschungszen- 
trum Karlsruhe and to answer it the construction of the new, 
powerful array which could be able to simultaneously measure 
as many shower characteristics as possible seems to be nec- 
essary . All particle components: electromagnetic, penetrating 
and hadronic carry the information needed to find the CR com- 
position on the event-by-event basis. The array KASCADE 
(KArlsruhe Shower and Core Array DEtector) started to 
produce physical results in late '90. Apart from experimental 
results there was another great achievement of the Karlsruhe 
group - the shower simulation program CORSIKA (COsmic 
Ray Simulations for KAscade) [6]. The importance of the 
CORSIKA code was not only in its detailedness, exactness 
and completeness, but mostly in that that it was made widely 
available, user-friendly and clearly documented. The helpful 
assistance, extensive debugging with continuous significant 
improvements of physics involved, made the CORSIKA a kind 



of a standard tool to be used by different experimental groups 
and theoreticians to explore and compare results on cosmic 
rays also in the energy regions far from initially imposed, 'the 
knee' region around 10^^-10^^ eV. These energies are not only 
of special interest for studies of the acceleration/confinement 
of CR in the Galaxy, but it also coincide, almost, with the 
highest accelerator energies available at that time. The high- 
energy interaction models which are essential for the EAS 
development were tested and adjusted to the SPS and Tevatron 
data. 

Originally CORSIKA consists of only one high-energy 
interaction model, the Dual Parton picture inspired model of 
Capdevielle, HDPM [7]. Then other models available on the 
marked were build-in to the CORSIKA structure as an options 
to be chosen by the user, simply by switching an appropriate 
flag in the steering cards. The history of major improvements 
is listed below. 

1989 CORSIKA 1.0 

SH2C-60-K-OSL-E-SPEC (Grieder, 1980) 
main structure, isobar model 
ESKAR (HDPM) (CapdevieUe, 1987) 

high-energy interaction 
EGS4 (Nelson et al. 1985) 

electro-gamma shower 
NKG (CapdevieUe, 1989) 

analytic EM subshowers 

1994 CORS1KA4.006 

GHEISHA (Feselfeld, 1985) 
VENUS (Werner, 1993) 
1997 CORSIKA 5.20 

SIBYLL (Fletcher, Geisser et al. 1994) 
QGSJET (Kalmykov et al. 1993) 
DPMJET (Ranft, 1995) 
2000 CORSIKA 6.00 

NEXUS (Drescher et al. 2001) 
UrQMD (Bleicher et al. 1999) 
2004 CORSIKA 6.20 

FLUKA (Fasso, Ferrari et al. 2001) 
2006 CORSIKA 6.535 

EPOS (Werner, Liu and Pierog 2006) 
(detailed references can be found in Ref.[8]) 

The most actual (June 27, 2008) version CORSIKA has a 
number 6.735 and consists of: 

EPOS - 161 Fluka 2006 

DPMJET - 11.55 GHEISHA^003d 
NEXUS - 397 URQMD - 1.3cors 

QGSJETOIC HERWIG 6.510 

QGSJET - n - 3 
SIBYLL 
VENUS 4.12/5 
The overlapping the time periods of the weakening of the 
CR and HEP relations and the proliferation of the CORSIKA, 
could be a pure coincidence, but it do not have to. In my, 
very personal opinion, the development of the CORSIKA, 
paradoxically, slowdown, and even, to some extent, suspend 
further CR driven development of high-energy interaction 



modeling. The multiparticle production codes introduced in 
CORSIKA overwhelmed all the others. Consecutive progress 
(see, e.g., [9]) in these particular codes made by the original 
Authors is important, but the last years show that the solution 
is still far-off. The solution means the correct description 
of EAS with the lack of contradictions from the accelerator 
measured characteristics. For recent example see, e.g., the 
paper by the KSCADE CoUaboration Ref. [10]. 

III. The simple analytic solution 

In this paper we would like to study EAS initiated by the 
particles (protons) of energies of 10^* eV and above. Such 
events are rare, thus the surface detectors to study them has 
to be spread over large area (if one wants to use the surface 
detector, there are other techniques available). The situation 
is quaUtatively different then the one around 'the knee' where 
the KASCADE measured each shower in details for its hadron, 
muon, and electromagnetic constituents. Such multicomponent 
and accurate measurements could give a profitable information 
about the interaction properties, but, on the other hand, this 
information is sometimes hard to handle, as it was just men- 
tioned. For very high energies only some characteristics of the 
shower can be measured with the enough accuracy. The most 
important is the shower longitudinal development, (given ex- 
perimentally as, e.g., the distribution along the shower axis of 
the fluorescent light emitted when the charged shower particles 
excite air molecules). This distribution could be measured to 
some extent only. The normahzation (total number of particles) 
which is related closely to the primary particle energy, and first 
moment (given usually as Xmax - the position of the maximum 
number of the particles in the cascade) are the parameters 
available for further study. The recent measurements of PAO 
[11] and HiRes [5] together with the older from Fly's Eye [12] 
and Yakutsk [13] gives the whole information of the average 
value of .Tniax as a function of estimated primary particle 
energy one can get. 

The interpretation of shower data in 'the knee' region needs 
simulation programs as complicated as CORSIKA to utilize 
the gathered information fully. To explore only the average 
aJmax data we can use a much simpler code. 

If we denote by iV^p^^^j {E, x) number of particles of the 
type p, -K,K ox jjL crossing the depth x with energies between 
E and E + dE, then N should evolve traversing the slab of 
matter of the depth dx 

- decreasing according to some probability of interact or 
decay within dx, 

- and possibly increasing by the average amount of parti- 
cles produced by the higher energy particles entering the 
slab. 

The only parameters of the respective system of integral- 
differential transport equations are probability of interaction 
(decay) and the inclusive energetic characteristics of multipar- 
ticle production processes. The decay constants and branching 
ratios needed for our purposes are known very weU so we'U 
concentrate on interactions hereafter. 



The interaction probabilities are given by the cross sec- 
tions (we'll also leave aside the problem of geometry and 
atmosphere modelling). They are of course not measured for 
the energies of our present interests, so we must extrapolate 
the low energy values. Changing the cross sections we can 
obviously control the rate of the shower development and 
move a shower maximum where it we wish it to be. But there 
are some constrains, and in recent years they become quite 
strong. We will discussed this point later. 

The inclusive energy distributions are known up to SPS (and 
Tevatron) energies and the knowledge is not as good as we 
wish it to be. Lower energy experiments (around \/2 « 20 
GeV) with stationary target produced more precise data and, 
what is more important, they cover the very forward region, 
which in fact controls a development of a cascade in thick 
media. The information we get is limited and we need to 
follow some more or less elaborate models to extrapolate it. 

A. HDPM 

The high-energy interaction model used in the CORSIKA 
program from the very beginning was the so-called HDPM 
model [7]. This phenomenological parametrization of the 
available data in the form expected by the two-chain structure 
expected in the Dual Parton jet hadronization pictures. In 
the proton-proton case it assumes on average that secondary 
particle rapidities are distributed according to two chains (jets) 
described by Gaussian (in rapidity). The widths, heights and 
positions of these Gaussians, as shown in the Fig. [1] are three 
most important parameters of the model. They are adjusted to 
the data up to SPS energies as shown in Fig.|2]and extrapolated 
smoothly. The parameters describe roughly the inelasticity and 
the average multiplicity. As the energy increase both increase 
as it is shown in Figs. 3 and 4. 




cms 



Fig. 1. Definition of three main parameters of HDPM. Tlie dashed (blue) 
components shows the contiibution introduced by the air target nucleus 
additional chains. 




Fig. 2. The SPS and Tevatron data on (pseudo)rapidity distribution compared 
with HDPM prediction with standard values of the model parameters 



B. Wdowczyk and Wolfendale scaling breaking model 

Another interaction model which we wont to examine in the 
present work is the scaling breaking model of Wdowczyk and 
Wolfendale (WW) ([14]). It has been proposed to described the 
CR data at the beginning of '70. Formally it is a generalization 
of the Feynman scaling idea given by: 
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where ^/i is the interaction c.m.s. energy E, p|| and p_L 
are energy, longitudinal and transverse momentum of created 
particle and Pmax is the maximum momentum which can be 
taken by the particle of the particular type (mass). 

The scaling, as it is widely known, was suggested in Ref.[3f] 
The WW modification is in the additional term 



2E 



f 



P\\ 



So 



,Pa 



(2) 



dp\\dp±_ \SoJ VPmax 

If the parameter of the model a is equal then the Feynman 
scaling is restored. For a = 0.25 the interaction multiplicity 
increases as E^^^ and, in the sense, this value can be treated 
as the upper, thermodynamical limit. The value of a originally 
introduced to the WW model is 0.13. This value is based on 
the interpolation of the xp ~ p\\/pmax distributions between 
sqrt{s) « 10 GeV and ISR energies. The increase of the 
central rapidity density reported in Ref.[15] suggests a = 
0.105. There are evidences from the high energy cosmic ray 
data that alpha could be even as big as 0.18. 



'It is interesting to note that the same idea appeared in the CR 
physics twenty years earlier in the paper by Heitler and Janossy: 
"#(E', E) = *(^) ^ shaU be a function of only." (where 
E is the incoming particle energy, E energy of the secondary particle and <& 
is the probability density of producing particle E). 



For the EAS description the WW model in its version of 
the mid '80 was improved by introducing partial inelasticities 
k (s, So) the slowly changing functions, which can give a 
better description of the production of different kind of sec- 
ondaries. The model predicts, for example, the increasing role 
of the production of the barions, and this was reaUzed by the 
power-law correction factor k with an index of 0.042. 
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The agreement of the WW model description of the SPS 
data was shown e.g. in Ref. [15]. 

IV. Results of HDPM and WW models at high 

ENERGIES 

Results of the comparison of the standard HDPM an the 

WW model at very high energies is given below. For the value 
of a we used for the moment 0.105 as suggested in Ref. [15]. 

The change of the inclusive distributions, rapidity for 
HDPM and x p for the WW model, as the laboratory energy of 
projectile proton increases from 10^'* to 10^^ eV is presented 
in Fig. 3. 





Fig. 4. Average inelasticity in p-Air interaction as a function of the proton 
energy for HDPM (left) and WW model (right). 

The differences in interaction inelasticity and in average 
multiphcities are shown in Figs. 4 and 5. It is worth to 
noticed that they are not very big. It is hard to speculate 
on the basis of these figures about the resulting difference 
in Xmax position. Surprisingly it is substantial. It comes from 
the behaviour of the inclusive x (or y and correlation of it 
with p±) distribution if the forward region rather than from 
global interaction characteristics. The WW model even with a 



= 0.105 brakes the Feynman scaling stronger than the HDPM 
change of the respective Gaussian widths. 

With the help of fast analytic program the average position 
of maximum of the shower can be calculated very fast. Results 
are shown in Fig. 6 for both models. The thin dashed lines 
shows the results of other (SIBYLL, DPMJET, QGSJET) 
CORSIKA models (the same is shown by the dashed lines 
in Fig. 5). These CORSIKA results in fact are well known 
from literature. 

The lack of the direct data which could help to solve the 
problem about the real nature of the forward fragmentation. 
We can try in this paper to use the cosmic ray data and the 
assumption about the 'pure proton' composition above the 
ankle to look for a possible solution. 

We will try to change an extrapolation of the model param- 
eters logarithmically as, e.g., in Ref. [16]. 

f 1 E< IPeV 

fne^{E) = Md{E) I ^ ^ ^^^^ _ M|2Z^ ^ ^ ^p^^ 

(4) 

This mean that, up to the energies of SPS (roughly). 



where the models have been adjusted to the data, we do not 
change anything and for higher energies the approximations 
are corrected by a slowly varying factor determined by the 
value of the correction Siq almost at the end of CR data, at 



IQi^ eV. 

This correction can be applied to multiparticle production 
model parameters as well as to the cross section values. But 
this has to be done with care. 




Fig. 5. Average multiplicity in p-Air interaction as a function of the proton energy for HDPM (left) and WW model (right). 
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Fig. 6. The prediction the proton induced shower Xmax depths calculated with HDPM (left and WW model (right) in a comparison with different measurement 
data. 



V. Recent progress in cross section description A. Hadron-pwton scattering cross section 

Recent years brought the significant progress in the theoret- 
ical description of the inelastic (elastic and total) cross section. 
The cross section involved in the EAS development is of It is based on the optical picture 
course hadron-nucleus cross section. Before we discuss this 
complex case we would like to look closely at the interaction 
with single proton. 
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The phase shift x is related to the scattering amplitude by the 
two-dimensional Fourier transform 
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Using the optical analogy one can interpret the 1 — e'-^^'^-' 
function as a transmission coefficient for a given impact 
parameter. Considering two colliding object we can assume 
(for pure absorptive potential) 



X{b) = iuj{b) 



Kab J d^h' paih)pb{h + b'), (7) 



with 



W{b;p) = ^{pbfKsifib) 



(9) 



give very good description of pp and pp data. Assuming the 
vector meson dominance and the additive quark model, it 
could be used with the same parameters also for photon-proton 
and photo-photon scattering cross section calculations. 
Another parametrization 

X(s,6) = (A(s)+^)a;(6,s) (10) 

was proposed by Perez-Peraza and collaborators in Ref. [18]. 

u}{b, s)= C {El Ko{ab) + E2 Ko{ab) + E3 K,i{ab) + 

Ei K,,{ab) + b [E5 Ki{ab) + Ee Ki{(3b)]} (11) 

was fitted to differential elastic scattering data. The consis- 
tency of ISR, SPS and Tevatror0 cross section within the 
framework of adopted model makes the eventual confidence 
band narrow. 

In the paper by Ishida and Igi [19] the cross section of K^p, 
TT^p, pp and pp the scattering amplitude was parametrised in 
the form 
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^Tevatron cross section is assumed to be of 74.21±1.35, tlie weighted 
aritlimetic mean of tlie E710 (72.8±3.1 mb), CDF (80.3±2.3 mb) and a E811 
(71.7±2.0 mb) values 



where ph is a particle's "opaqueness" (the matter density 
integrated along the collision axis). 

In the serial of papers by Block and coworkers [17] the 
approximation of x of the form inspired by QCD 

X{b,s) = £,qq{s,b)+^qg{s,b)+£,gg{s,b) = 

= i [aqqW{b] p,qq) + CTqgWib] jJiqg) + 'JggW{b\ fJ.gg)] (8) 



which, for high energies, leads to the saturation of the Froissart 
bound 

a ~ Blog2(s/so) = (J^^C2\og'{s/so) (13) 

The universality of the value of B found in Ref. [19] 
0.289±0.023, 0.351±0.036, 0.37±0.21 for pp (pp), 7r±, and 
K^, respectively, gives additional evidence that proposed 
picture is correct and extrapolations are thus strongly justified. 

The very similar and self-consistent description of data on 
charged pion and proton (antiproton) projectiles has been given 
by Block and Halzen in Ref [20]. 

cr(s) = CO + Ci log ( — ) + C2 log^ ( — ) + /3p' ( — ) 

Vm/ Vto/ V771/ 

(14) 

With some new tool of calculating 'the best fit' the necessity 
of log^ component is confirmed. 

The log'^ character of the cross section rise was argued for 
by COMPETE Collaboration in Ref. [21]. 

We would like to remind here briefly the result of Ref. [22] 
from 1998. 

The cross sections for K^p, n^p, pp and pp interactions 
were parametrized assumed geometrical scaling 
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Fig. 7. Proton-proton total cross section calculated using parametrisation 
Eq.jTs) compared with the data from accelerator measurements an recalcu- 
lated cosmic ray attenuation data points. Dashed and dotted lines represent 
the inelastic and elastic contributions. Dotted and dash-dotted (coloured) lines 
are the results from other well known approximations [23], [24], [25]. 
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where omega was calculated as a convolution of colliding 
hadron matter distributions described with the help of only one 
parameter {mTi,mK, rup) for each interacting particle type. 
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Values of A and CTinci are found: 

0.077 ln(s/so) 



Xis) 



1 + 0.18 ln(s/so) + .015 In^is/so) ' 



CTi„ci(s) = 32.4 - 1.21n(s) + 0.211n2(s) 
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(sq = 500 GeV^). The resulting cross sections are given in 
Fig |7] together with other approximations by Block and Cahn 
[23], Durand and Pi [24] and Honda parametrization using 
Akeno data [25]. 

B. Proton-nucleus and nucleus-nucleus interaction cross sec- 
tion 

The scattering of particle on the close many-particle system 
(nucleus) can be treated as a superposition of individual 
interactions each with a specific phase shift. The overall phase 
shift for incoming wave is a sum of all the two-particle phase 
shifts. 



(19) 



XA{b, {d}) = - d,) 

Eq. (fTgj l is the essence of the Ref.[26] and defines the Glauber 
approximation. 
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Fig. 8. The relation between pp and P-Air cross sections calculated with 
exact Glauber formula and other dependencies used to convert cosmic ray 
EAS attenuation data to pp cross section. 

On the other hand, the scattering process can be treated as 
the single collision process with its own nuclear phase shift 
Xopt(^) To get the consistency with Eq.(fT9l) it is required 



piXopt(ft) _ 



m{d}re 



To go further with the calculations of Xopt a commonly used 
assumption has to be made. If we assume that the number 
of scattering centers (A) is large and the transparency of the 
nucleus as a whole remains constant then 



e^x(^{d})\) (20) 

what defines the relation of the individual projectile-nucleon 
and overall projectile-nucleus oppacities. 
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where pA is the distribution of scattering center (nucleon) 
positions in the nucleus (^ gj). 
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Fig. 9. Proton-Air interaction cross section calculated according to exact 
Glauber formula and the values used by DPMJET QGSJET and SIBYLL 
models of CORSIKA program. 

And finally, assuming that the individual nucleon opacity 
|1 — e*-^'-'''! is a very sharply peaked compared with pA then 
with the help of the optical theorem the simple formula can 
be found 
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where the last equality holds in the large A limit (certainly 
Eq.(l22li cannot be used for A — 1) This result is often but not 
quite correctly called "the Glauber approximation". As it has 
been shown, the original Glauber assumption given by Eq.(fT9l) 
is here supported by small nucleon size and a large value of 
A. 

As we saw for very high energies nucleons are quite big 
objects and it is expected that at least the last approximation 
could be questioned. 

We have performed respective calculations. The results 
is shown in the Fig|8] as a relation between pp and p-Air 
interaction cross sections. It is shown there with the other 
used in the literature. The importance of this relation is that it 
allows one to get the pp cross section from the cosmic ray data 
on EAS attenuation length measured experimentally which is 
related to p-Air interaction cross section. 



TABLE I 

Predicted values of pp cross section at LHC energies. 



Author 


Year 


Ref. 




Honda 


1993 


[25] 


110.4 


Wibig and Sobczyiiska 


1997 


[22] 


102.5 


Cudell et al. 


2002 


[21] 


111.5±1.2 


Velasco et al. 


1999 


[29] 


104.17±4.4 


Perez-Peraza et al. 
Block et al. 


2005 
1999 


[18] 
[17] 


108.27 
108±3.4 


Block and Halzen 


2006 


[20] 


107.3±1.2 


Ishida and Igi 


2007 


[19] 


109.5±2.8 



The consistency of our cross section description is shown in 
Fig. m as the p-Air cross section calculated with exact Glauber 
formula Eq.(fT9]l and pp inelastic cross section parametrization 
of Eqs.dTTl [TSl l. There are also shown cross sections adopted 
by various high energy interaction models in CORSIKA. 

Concluding the discussion on cross sections, we show in 
the Table H] some predictions of various authors concerning he 
pp total cross section for the LHC energy ^/s = 14 TeV. 

As it is seen, all the estimated values are very close and, 
with high degree of confidence, it can be concluded that the pp 
cross section predicted for the LHC energies is expected to be 
equal about 108 mb (within few millibarns 'error box'). The 
result 102.5 mb is excluded by the most later fits. However, 
it should be underline here again that this result was obtained 
by the exact Glauber formulas while the 108 mb prediction is 
influenced by the HAS data obtained using the transition from 
p-Air to pp made with the help of point nucleon approximation 
or even with multiple scattering approach [30]. 

VI. 'Protons only' and HDPM results on Xmax 

As if is shown in Fig. 6 (left) the pure proton flux around 
the ankle is excluded when HDPM with default parameters 
is used. The analytic shower development allows us to test if 
the correction factor of the form of Eq.© applied to a model 
parameter can make the 'pure proton' hypothesis acceptable. 



We can try first to change the interaction cross section. It 
is possible to get some nice result as it is seen in Fig. 10. 

But this results needs cross section as the one shown in the 
small inserted plot in Fig. 10 what is certainly unacceptable. 

Applying the correction given by Eq.(|4|l (with Sig not 
greater than ^ 3) to the width of the Gaussian (Fig. [T]l, the 
average multiplicities (taking care not to exceed the available 
energy, inelasticity coefficient can not be bigger than 1), or to 
average p± the 'pure proton' Xmax can not get close to the 
measured values. 

850 I ^ 1 




E (eV) 

Fig. 10. HDPM 'best' results for 'pure protons' with the cross section 
increased. The inserted plot shows a respective increase of a^^^ in the very 
high energy region. 

Concluding: there is no way to adjust the HDPM parameters 
in the way that the position of the shower maximum agree with 
measurements for pure proton composition at and above the 
ankle. 
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Fig. 1 1 . Position of the shower maximum as a function of primary proton energy calculated with the WW model for constant parameter a=0 (left) and 0=0.07 
with slightly increased multiplicity and other parameters (right) compared with data. Small inserted plots shown the average multiplicities respectively. 



VII. Xmax RESULTS FOR 'PROTONS ONLY' WITH WW 
MODEL OF HIGH ENERGY INTERACTION 

The situation for the Wdowczyk and Wolfendale model 
is different. The proposed strong violation of the Feynman 
scaling came as a result of analysing the variety of cosmic ray 
data rather than some suggestions from simplified theoretical 
picture of jet fragmentation where the theory (QCD-like) has 
still problems with calculating non-perturbative effects. Re- 
sults shown in Fig. 6 give a hope that the WW parametrization 
could be able to describe the position of the shower maximum 
at very high energies with pure proton primary spectrum. 
The important question here is if the eventually adjusted 
parameters will give an acceptable characteristics for the single 
pp interaction. 

With the help of the fast analytic program the proposed 
changes to the average p^, multiplicity and inelasticity has 
been tested, but the main interest was put to the value of the 
WW model parameter a. Fig. 6 shows results for q;=0.105. In 
Fig. 11 results for a = (Feynman scaling) and a = 0.07. 

Values of .Tmax for Feynman scaling are much to deep in the 
atmosphere, and looking into detail of the average multiplicity 
of produced secondaries (shown in Fig. 11 as inserted small 
plots) one can find that they are far too small in comparison 
with any expectations. 
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Fig. 12. The 'best' Xmax result of Wdowczyk and Wolfendale model for 
'pure protons'. The inserted plot shows the value of the model parameter a. 



provement can be obtained with a changes of the average 
multiplicities according to the form of Eq.(|4|. However, the 
multiplicities can not be upraised much without a change of 
secondary particles energy distribution because the limited 
energy available. The inelasticity calculated by dividing the 
integrated secondaries energies by the incoming energy can 
reach at most 1 (the elasticity defined as the fraction energy 
carried of the most energetic particle is of course smaller 
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Fig. 13. Average inelasticity (left) and multiplicity (right) of the WW model 
eventually adjusted to a'max data in Fig |12l 



Much better results can be achieved allowing the slow 
increase of the alpha from 0.07 at lO^^eV to 0.134 at lO^^eV 

Average inelasticity and multiplicities at very high energies 
of such 'best fit' of pp interaction scaling breaking are shown 
in Fig. 13. 

The rise of a can be translated to the behaviour of the 
effective WW scaling factor (s/sq)" modifying the inclusive 
xp distributions. This is shown in FigdU It eventually goes 
like the factor of 0.13 found in original Wdowczyk and 
Wolfendale fits from more than a quarter of century ago. 
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Fig. 14. The WW-scaling factor (s/sq)" of the lit shown in Fig. 1121 

VIII. Conclusions 

We have analyzed the possibility of the pure proton com- 
position of ultra high energy CR proposed recently by the 
PAO and HiRes experiment. This assumption contradicts the 
conventional wisdom on the rate of the shower development 



measured by different big apparatus for some time and known 
as the depth of the shower maximum energy dependence. The 
interaction models incorporated in the EAS simulation code 
CORSIKA all give the composition enriched by heavies with 
average logarithm of the mass in about the middle between 
proton and iron. It is clear that the pure proton composition 
needs a change of the interaction model to give the self- 
consistent picture of the nature of the CR flux. The way of 
the changes has been shown. There is necessary to violate the 
Feynman scaling very strongly. We have found that the scaling 
breaking model of Wdowczyk and Wolfendale is well suited 
to study this. The WW model parameter a which describes 
the Xmax PAO data with the pure proton CR is close to the 
value found about a quarter of century ago by Wdowczyk 
and Wolfendale in the original papers. The spread of data 
points and reported error bars with possible systematics do 
not allow us to perform fits much more precise that the 
one presented in this paper It can be said that the total 
multiplicities could be slightly below expectations, but this 
is related mostly to the central rapidity region and the WW 
model is not supposed to specify just this region exactly. The 
average inelasticity could be assumed constant, but even slight 
decrease is possible. We have shown also that there is no way 
to find out the modification of Dual Parton inspired models 
(like HDPM) to be adjusted to the pure proton flux and Xmax 
as measured for giant EAS. If one wants to have protons 
only than for the interpretation of the shower data, energy 
calibration etc. another interaction model has to be introduced 
to the CORSIKA repository. 
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